Atom laser divergence 
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We measure the angular divergence of a quasi-continuous, rf-outcoupled, free-falling atom laser as 
a function of the outcoupling frequency. The data is compared to a Gaussian-beam model of laser 
propagation that generalizes the standard formalism of photonic lasers. Our treatment includes 
diffraction, magnetic lensing, and interaction between the atom laser and the condensate. We find 
that the dominant source of divergence is the condensate-laser interaction. 

PACS numbers: PACS numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj, 04.30.Nk 
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A dilute gas of atoms condensed in a single quan- 
tum state of a magnetic trap [[j} |[ ||, |], ||, || is the 
matter-wave analog of photons stored in an optical cav- 
ity. A further atomic parallel to photonics is the "atom 
laser" [0, ||, ^|, 10 , a coherent extraction of atoms from 
a Bosc-Einstcin condensate. Atom lasers are of basic 
interest as a probe of condensate properties flll| and 
as a coherent source of atoms [ fL2| . Furthermore, just 
as optical lasers greatly exceeded previous light sources 
in brightness, atom lasers will be useful in experiments 
that simultaneously require monochromaticity, collima- 
tion, and intensity such as holographic atom lithogra- 
phy (l3| , continuous atomic clocks H] , and coupling into 
atomic waveguides jjjl . 

Initial work on atom lasers has demonstrated several 
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types of coherent, pulsed output couplers [J 
as a coherent narrow-band coupler jlfl, 06 
type employed in our work. In this Letter, we address 
the nature of propagation of an atom laser outcouplcd 
from a condensate. Our experimental data show that 
the laser beam is well characterized by a divergence an- 
gle. We measure this angle versus radio-frequency (rf) 
outcoupler frequency, which chooses the vertical extrac- 
tion point of the atom laser from the condensate (see 
Fig. |l|a). In choosing the extraction point, one chooses 
the thickness of the condensate to be crossed by the ex- 
tracted atoms, as well as the width of the atom laser 
beam at the extraction plane. To interpret this data, 
we use a formalism that generalizes the ABCD matri- 
ces treatment of photonic lasers [J17], [18[ . This treatment 
allows us to calculate analytically the divergence of the 
laser due to diffraction, magnetic lensing, and interac- 
tions with the condensate. We find that, for our typical 
experimental conditions, the divergence of the laser is pri- 
marily due to interaction between the atoms in the laser 
and the atoms remaining in the condensate. We describe 
this interaction as a thin lens. Note that the existence 
of such an interactive lensing effect is in stark contrast 
to a photonic laser, since photons do not interact. In- 
teractions were estimated to be similarly important for 
pulsed atom lasers J?], [|. In our case, the divergence is 
also magnified by a thick-lens-like potential due to the 
quadratic Zccman effect of the magnetic trapping fields. 



The technique we use to obtain Bose-Einstein con- 
densates with 87 Rb is described in detail elsewhere [ p^9[ . 
Briefly, a Zeeman-slowed atomic beam loads a MOT in 
a glass cell. Typically 10 8 atoms are transferred to a 
magnetic trap, which is subsequently compressed to os- 
cillation frequencies of ui x — co z = 2n x 144 Hz and 
Lu y = 27r x 9 Hz= \uj z in the quadrupole and dipolc 
directions, respectively, where z is vertical. The Ioffc- 
Pritchard trap is created with an iron-core electromag- 
net, with a typical quadrupole gradient of 11.7 T-m , 
and an uncompensated dipole bias field of Bq =5.4 mT. 
A 40 s, rf-induced evaporative cooling ramp in the com- 
pressed trap results in condensates of typically 5 x 10 5 
atoms, with a 15% rms shot-to-shot variation. 

Atom lasers are created by applying a rf field at about 
38.6 MHz, to transfer condensate atoms from the trapped 
\F,itif) = 1 1 , — 1} state to the weakly anti-trapped 1 1 , 0) 
state, which also falls under gravity (see Fig. |l]a). The 
rf field is weak (approximately 0.1 //T) and applied for 
a relatively long duration (t oc = 10 ms) |16|. There is 
no significant coupling to the 1 1, — b-1) state because the 
F = 1 sublevel transitions are split by 0.8 MHz due to 
the nonlinear part of the Zeeman effect at 5.4 mT. To 
measure the spatial distribution of the atom laser, we 
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FIG. 1: (a) Schematic representation of outcoupling from the 
trap potential Vji _j\ to the untrapped Vn o) potential via an 
rf photon at cj r f . The vertical direction is z. The curves shown 
include magnetic, gravitational, and mean field potentials. 
(b) The wave-vector k b s for the absorptive imaging light 
beam at an angle cf> with respect to the weak axis y of the 
condensate. The horizontal coordinate in the imaging plane 
is r. In both images, the condensate is cross-hatched. 



take an absorptive image with a pulse of resonant light 
6 ms after the moment when the rf outcoupling ends and 
the trap is turned off. As depicted in Fig. |l]b, the image 
is taken at a <f> — 55° angle from the weakly confining y 
axis of the trap. Figure @a shows a typical image of an 
atom laser. 

The gravitational sag <?/w 2 shifts the entire conden- 
sate from the center of the magnetic trap to a region 
where iso-field surfaces are approximately planes of con- 
stant height zq across the condensate p}. The relation 
between rf frequency luq + S T { and coupling height zq is 
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where A = MgR z /H is the spectral half- width of the 
condensate, M is the mass of the atom, g is grav- 
itational acceleration, and R z is the Thomas-Fermi 
(TF) radius |2Q| of th e condensate along z. The ra- 
dius R z is y/2n/Mu}1, and fx is the chemical potential 
{hLu/2)ll5a 11 N) 2 / 5 (M[u/h)- 1 / 5 , where N is the number 
of atoms in the condensate, an = 5.67 nm is the s-wave 
scattering length between 87 Rb atoms in the |1, — 1) state 
pl[ , an d w 3 = id x LUyUj z . In Eq. (|y), we have chosen loq 
such that <5 r f = ±A when zq = =Fi? 2 . For our experimen- 
tal parameters, /z/ft <C A, so 5 T { is roughly linearly depen- 
dent on zq, with slope —Mg/h. Therefore, in choosing 
the coupling frequency, we choose the height within the 
condensate at which the laser is sourced. 

Short-term (ms-scale) stability of the bias field is ver- 
ified by the continuity of flux along the laser (Fig. pk). 
We find that the shot-to-shot stability and reproducibil- 
ity of the bias field, between each 80 s cycle of trapping, 
cooling, and condensation, is about ±0.4 ^T or ±3 kHz. 
This stability is sufficient to scan through the condensate 
spectral width of about 20 kHz. In one out of five runs 
(a run is a set of about 10 cycles), the data were not self- 
consistent, which we attribute to a larger (>l/xT) bias 
field fluctuation during that run. We could maintain this 
bias field stability, typically one part in 10 4 [f22| , by using 
either a low-noise power supply or a battery, since only 
1.3 A is used to energize the dipole coils during evapora- 
tion and outcoupling. 

We analyze the images (such as Fig. ya) to measure 
the flux and divergence of the outcoupled laser. Fig- 
ures ||b through gd shows the first step in the analy- 
sis, a series of fits to the transverse spatial profile of the 
laser, as measured at several heights z n , averaging across 
±140 /im. The column density profile is fit at each z n by 
p n {\ — r 2 /i? 2 ) 3 / 2 , where r is the horizontal coordinate 
in the image plane (see Fig. Up), p n is the peak column 
density, and R n is the width. This fit function would be 
the rigorously correct function for an atom laser without 
divergence [g3| or for a free condensate undergoing mean 
field expansion p4fl; here, we observe empirically that the 
fits are good. From the integral of the column densities 
at each height, we determine the output flux F with a 
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FIG. 2: Typical continuous atom laser output, (a) Ab- 
sorptive image of laser, after 10 ms of output coupling at 
38.557 MHz and 6 ms of free flight. The condensate (the dark- 
est area of image) is displaced from the beginning of the laser 
because of a magnetic kick separating atoms in the tuf = — 1 
and m,F = states, (b)-(d) Examples of absorption profiles 
of the atom laser, taken from the three regions boxed in (a). 
Each profile is fit to find a width R n and density p n . (e) The 
divergence angle of a single laser is found with a linear fit to 
the series of measured widths {/?„}. 



fit to F/ y^2g(z — zo), a form that assumes constant flux, 
purely gravitational acceleration, and a density that de- 
creases with the inverse of the classical velocity, valid 
when \z-z \ > £ g , where £ g = [h 2 /(2gM 2 )} 1 / 3 is the 
gravitational length. Finally, from the series {R n } and 
{z„}, we determine a geometric expansion angle by a 
linear fit (see Fig. ge). We will discuss below why one 
would expect a linear rather than parabolic shape for a 
laser falling under gravity. 

We repeated the above imaging, and analysis for atom 
lasers coupled at a variety of rf coupling frequencies. Fig- 
ure @ shows the averaged half-angle divergence versus 
outcoupler detuning, from 15 laser images. We see that 
the divergence clearly decreases at higher S r { , correspond- 
ing to lower initial heights zq. These divergence data will 
be analyzed in more detail below. The inset of Fig. H 
shows laser flux F versus <5 r f , and is compared to 
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where Fq is the peak flux, and Eq. (pi) defines the relation 
between zq and 6 T {. Equation (0) assumes the laser flux 
is simply proportional to the linear density of the con- 
densate at the coupling point. For the solid curve shown 
in the inset of Fig. H, we have used N = 4 x 10 5 to give 
A = 2tt x 9.1kHz and R z = 4.2/im. The peak flux is 
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FIG. 3: Divergence (half-angle) versus output coupler detun- 
ing from the condensate center. Experimental points are com- 
pared with the theoretical calculation (solid line) . The dashed 
line represents the same calculation, but excluding the effects 
of the condensate-laser interaction. Inset: Output flux (in ar- 
bitrary units) versus detuning, with the same frequency scale 
as the main figure. 



measured to be about 1 x 10 7 s _1 by the reduction in 
condensate number, as in llQ]. This is in agreement with 
the theory of (23], given the applied field strength and 
uncertainty in the atom number and field polarization 

There are several possible sources of divergence of the 
atom laser, including diffraction, magnetic lensing, and 
interactions both within the laser and between the laser 
and the condensate. In order to understand the diver- 
gence with a simple analytical model, we make several 
approximations: 1) that the interactions between atoms 
within the laser are not significant, valid in the low- 
flux limit; 2) that the roughly parabolic density pro- 
file of the atom laser can be approximated by a Gaus- 
sian; and 3) that we can use stationary solutions of the 
Schrodinger Equation with a paraxial-type of approxima- 
tion, in which the fast degrees of freedom [g3| are decou- 
pled from the slow evolution of the transverse degrees of 
freedom, as in (25). We follow a Gaussian optics treat- 
ment similar to that of photonic lasers [flTj : the spatial 
distribution in x and y described by the wavefunction 
$(x,y,t) = y Q exp[-iP(t) + ix 2 /2q x (t)+iy 2 /2q y (t)], 
where P(t) describes the overall phase and amplitude, 
and beam parameters l/q x . y (t) — i/w xy {t) + c XyV {t) de- 
scribe the widths w x and w y and curvatures c x and 
c v of the beam. The observable width is w(t) = 



(w 2 (£) cos 2 ((f>) + WySin 2 (0)) 1 / 2 , where 4> is the obser- 
vation angle (Fig. |l]b). The initial widths w 2 (0) = 
2{R 2 — 2q)/5 and w y (0) = w x (0)/X give the same rms 
width as would an initial Thomas-Fermi density profile. 
The beam parameters q x and q y follow an "ABCD 
law" similar to that for an photon laser beam: q' = 
(Aq + B)/(Cq + D), where the coefficients A, B, C, and 
D are the four elements in a matrix which transform ray 



vectors (x,p/h) to (x',p' /h) according to classical equa- 
tions of motion in the same potential. In the following 
paragraphs, we will calculate the ABCD matrices for in- 
teraction with the condensate, propagation with the mag- 
netic trap on, and free flight after the trap is turned off. 
Note that even though the ray matrices can be derived 
using classical equations of motion, their application in 
the Gaussian beam formalism includes diffraction. 

The mean-field interaction potential between an atom 
in the laser and the condensate is Ui(r) = g ip c (r), 
where 501 is the s-wave coupling strength 47rfi 2 aoi/M be- 
tween atoms in the |1,0) state and the trapped |1, — 1) 
state, and p c (r) is the condensate density. Here we use 
aoi = o-n- We calculate the action of this potential 
treating it as a lens, and using the thin lens approxi- 
mation that each trajectory is at its initial transverse 
position (xo,yo). In the Thomas-Fermi limit, the po- 
tential Ui(r) is quadratic, and thus gives an impulse 
(Ap x ,Ap y ) = (2fit 1 (x Q ,y Q )/R 2 ) (x ,X 2 y Q ) after an in- 
teraction time t\. The on-axis power of the lens is given 
for xq <C R x and yo <C R y , with which assumption we 
find t\ w 2(R Z + z )/g. This gives the thin lens ABCD 
matrix for the x direction 



Mi x (z a ) = jfj_ 



1 



J^^/2(R z + z ) 
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When applied to the beam parameter q x , the nontriv- 
ial term in Eq. (f|) is the wavefront curvature added to 
the beam. A similar ray matrix Mi v (zq) transforms q y , 
but with a curvature term multiplied by A 2 . For both 
Mi x (zq) and Mi y (zo), the curvature is positive for all 
zq, since the interaction is always repulsive. A positive 
curvature corresponds to an expanding wave, and thus 
the condensate with repulsive interactions (501 > 0) is 
always a diverging lens. 

When the atom laser falls [£(|, it evolves in the anti- 
trapping potential due to the quadratic Zeeman effect of 
the |1,0) state, C/qze(i") = ~-^bB 2 {t)/ -Bhf, where [Ib is 
the Bohr magneton, and £>hf = 0.4883 T is the hyperfine 
splitting in magnetic units. We can neglect y-dependent 
and higher-order terms, since they are several orders of 
magnitude smaller for the fields in our experiment, to 
get U QZE {r) « -ii B Bl/B UF - Mn 2 (x 2 + z 2 )/2, where 
£1 = 2-7T x 30.3 ±0.1 Hz. Below, we will return to evolu- 
tion in y. The classical motion of a particle in an inverted 
quadratic potential is given by hyperbolic functions. In 
the vertical direction, an elongation of the laser is evi- 
dent: we observe a length of 1.84±0.09mm, while with 
gravity alone, one would expect a length of 1.33 mm af- 
ter 10 ms of coupling and 6 ms of free flight. Including 
Uqze(y) and our measured trap parameters, we calcu- 
late 1.87 mm, in agreement with our observations. The 
ray matrix for the transverse x direction is 



M 2 {t 2 



cosh(f2t 2 ) jj^ sinh(S!£ 2 ) 
^sinh(m 2 ) cosh(nt 2 ) 



(4) 



where t 2 is the time of evolution, ranging between and 
t oc . This interaction is a thick lens, since f2i oc > 1, 
and thus there is sufficient time for the laser to change 
diameter and curvature during its interaction. During 
the same time, the beam parameter q y transforms by the 
free flight matrix 



M FF (t 2 ) = 



ht 2 /M 
f 



(5) 



measurements and model provides a basic tool for future 
experiments with atom lasers. 
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which is the fi — » limit of M 2 (t 2 ). Note that it is due to 
the acceleration in both the z and x directions that make 
the atom lasers well characterized by an asymptotic ex- 
pansion angle (see Fig. @e) : if there were no acceleration 
in x, the laser would have parabolic borders. 

The third and final transformation of the laser is free- 
flight expansion between turning the trap off and observ- 
ing the laser, described by M FF (t F ), where t F is the time 
of flight. To find the width of the laser at any position, we 
evolve the initial beam parameter q x (0) = —iw x (0) 2 and 
q y (0) = —iw y (0) 2 with the elements of the matrix prod- 
uct M FF (t F )M 2 (t 2 )M lx {z ) and M FF {t 2 +t F )M ly {z Q ), 
respectively. The angle is the ratio of the difference be- 
tween the observed rms beam size at t 2 = and t 2 = t oc 
to the length of the laser. Note that this theory has 
no adjustable parameters since the atom number, output 
flux, trap frequencies, interaction strength, and bias field 
have all been measured. 

Figure shows the calculated angle 9 in comparison 
with the data. The primary feature of this curve is its 
monotonous decrease with increasing 5 r { throughout the 
range of the data. This trend is due to the condensate 
lens, as is demonstrated by comparison to 9 without the 
transformation Mi (dashed line in Fig. 0). Simply put, 
a laser sourced at lower heights (greater <5 r f) interacts 
for less time with the condensate. The quadratic Zee- 
man potential J/qze acts to magnify the divergence of 
the condensate lens, increasing the slope in this region 
by roughly a factor of four. Out of the range of our 
data, there are two more salient features: (1) The angle 
9 decreases at (5 r f < —5 kHz. This is due to a reduction 
in initial width of lasers sourced from the very top of 
the condensate, since divergence angle is proportional to 
beam radius for a constant focal length. (2) The angle 9 
increases for S T { > 5 kHz due to diffraction. A combina- 
tion of diffraction and interactions imposes the minimum 
divergence possible on the atom laser, in our case approx- 
imately 6 mrad. 

In conclusion, we have measured the divergence of an 
atom laser. We demonstrate that in our case, interactions 
are a critical contributor to the observed divergence. The 
strong parallel between atom and photon lasers beams, 
both fully coherent, propagating waves, is emphasized 
by the success of a model obtained by generalization of 
the standard treatment of optical laser beams. The un- 
derstanding of atom laser propagation provided by our 
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